Rectification by imprinted phase in a Josephson junction 
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A Josephson phase shift can be induced in a Josephson junction by a strategically nearby pinned 
Abrikosov vortex (AV). For an asymmetric distribution of imprinted phase along the junction (con- 
trolled by the position of the AV) such a simple system is capable of rectification of ac current in a 
broad and tunable frequency range. The resulting rectified voltage is a consequence of the directed 
motion of a Josephson antivortex which forms a pair with the AV when at local equilibrium. The 
proposed realization of the ratchet potential by imprinted phase is more efficient than the asymmet- 
ric geometry of the junction itself, it is easily realizable experimentally, and provides rectification 
even in the absence of applied magnetic field. 

PACS numbers: 74.50.+r,74.45.+c, 74.78.Na, 74.25.Fy, 74.20.De 



Starting from the discovery of biological molecular mo- 
tors [l|, the ratchet effect has been demonstrated in 
many different physical systems, where rectification in 
the presence of external random or periodic forces with 
zero time average is induced bymeans of spatial or tem- 
poral asymmetries (see Ref. [2| for a review). Among 
other solid-state ratchet systems, superconducting ratch- 
ets have been realized - based on Abrikosov vortices 
01 • Josephson vortex ratchets were also studied, in long 
Josephson junctions (JJs) 0-0 and in specially engi- 
neered J J arrays [7]. Voltage rectification based on the 
Josephson phase change has been demonstrated in asym- 
metric [8, 9] and three-junction SQUIDs ^q|, [HI, and in 
annular JJs pj| with the asymmetric potential created 
by junction design, by inhomogeneous magnetic field [HI 
or by extra current biasing [6] . Josephson ratchets based 
on asymmetry of the drive rather than the potential it- 
self have also been realized [l^ . The nonlinear signal 
mixing of two driving forces was also shown to be ca- 
pable to control transport in different deterministic and 
Brownian ratchet devices 
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In this Letter, we propose a ratchet based on an inho- 
mogeneous phase change along the planar JJ. The sim- 
plest practical realization of such phase distribution can 
be realized by pinning an Abrikosov vortex (AV) nearby 
the junction (see Fig. 1). An AV can be inserted into the 
sample, e.g., by field cooling |17|], or by passing a large 
bias current through the system jl8| . Once in the system, 
the location of the AV can be controlled by e.g. appro- 
priately directed transport current [l^- However, such 
current would affect the J J as well. A more elegant way 
of nucleating, as well as manipulating the AV is through 
the use of an electron beam, demonstrated experimen- 
tally by Ustinov et al. [19]. Here, we keep the AV from 
penetrating into the junction, since its magnetic field is 
known to strongly alter the properties of the junction 

Our idea is fairly simple. In the presence of an ap- 
plied ac current, an off-center location of the AV creates 
an asymmetric phase imprint on the junction and, con- 



sequently, an asymmetric potential for the motion of a 
Josephson fluxon along the junction. This in turn gen- 
erates a net dc voltage, i.e. rectification, in a broader 
current and frequency range than earlier ratchets based 
on geometric asymmetry of the junction. Moreover, a 
clear advance of the proposed ratchet is that it operates 
in the absence of applied magnetic field, whereas most 
known superconducting ratchets are not. 

We consider a thin superconducting stripe with a nar- 
row metallic junction (see Fig. [T]). The location of the 
AV is predetermined by a hole which acts as a pinning 
site. To take into account the Josephson coupling and 
the phase field of the AV inside the junction, we use 
the Ginzburg-Landau (GL) formalism with a Lawrence- 
Doniach extension [2l|]. The resulting, modified time- 
dependent GL equation [22] takes the form 
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coupled with the equation for the electrostatic poten- 




FIG. 1: (color online) An oblique view of the system: a super- 
conducting stripe (width length L and thickness d <^ X,^) 
with a central Josephson junction (width dj). An Abrikosov 
vortex is trapped in a hole (of size a) located at a distance 
dx from the junction (^-axis) and dy from the x-axis. The 
current is applied via normal-metal contacts (labeled 1 and 
2) and the output voltage is measured at a distance 6 away 
from these leads. 
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FIG. 2: (color online) (a) I-V characteristics of the sample 
with L = 4 /im, w = 2 /im and dj — 20 nm without (solid 
curve) and with the Abrikosov vortex (symbols) at a distance 
dx away from the junction (dy — 0). (b) The critical current 
jc vs. for two sizes of the sample (symbols) together with 
the analytical result (solid line, see text). 
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FIG. 3: (color online) Voltage vs. time response of the sample 
at = 0.03jo for dx = 500 nm, upon transition to the 
resistive state (see Fig. EJja)). Panels (1-4) show snapshots of 
the phase (left) and density (right, just the junction area) of 
the order parameter at the times indicated in the V{t) curve, 
and illustrate the motion of the Josephson vortex along the 
junction. 



tial A(f = div(^[?/^*(V — iA)ip])^ where A is the mag- 
netic vector potential. These coupled non-linear differ- 
ential equations are solved self-consistently using Euler 
and multi-grid iterative procedures. The last term in Eq. 
(1) (with A = /^^^/2 ^xdx and /i the ratio of the mass 
of the Cooper-pairs in the superconducting and normal- 
metal regions) describes the Josephson coupling across 
the junction of width dj, and exists in the calculation 
just for the points bordering the junction on either side 
(indexed ±(ij/2). Here, the order parameter is scaled 
to its value at zero magnetic field i/jq, distances to co- 
herence length ^(0), time to tgl = Trh/SkBTcU^ vector 
potential to cfi/2e^(0), and the electrostatic potential to 
(fo = h/2eTGL- The material parameter F = 2re'0o/^ 
(with Tg being the inelastic scattering time) and u are 
chosen to be 10 and 5.79, respectively [22|. In the sim- 
ulations we take ^(0)=10 nm and A(0) = 200 nm, which 
are typical for Nb thin films [23|. Using the normal- 
state resistivity p = 18.7 pQcm for such films we obtain 
TGL ^ 2.69 ps and (po = 0.12 mV at T = 0.9Tc which is 
the considered temperature in our simulations. Neumann 
boundary condition is applied at all sample boundaries 
(including the hole), except at the current contacts where 
we use = and \/(f\n = — j, with j being the applied 
current density in units of jo = ch/ IQeTr^X^^. 

We begin our analysis by demonstrating the properties 
of a Josephson junction {dj = 20 nm) with symmetrically 
imprinted phase shift along the junction (i.e., dy = 0). 
Figure [Sfa) shows the time-averaged voltage vs. applied 
current (I-V) characteristics of the sample with length 
L = 4 jam and width w = 2 pm without (solid black 
curve) and with the AV (symbols) at different distance 



dx from the junction. The I-V curves exhibit a strong 
dependence on the vortex position - the critical current 
jc for the transition to the resistive state decreases as 
the AV is placed closer to the junction and the imprinted 
phase difference = 2arctan(tL'/2(ix)/7r increases. Fig. 
[2](b) shows the calculated jc as a function of A^ (open 
dots) together with the known analytic expression for 
the critical current of the junction for given phase differ- 
ence on the junction jmax/jc = sin(A^/2)/(A^/2) (solid 
line). The latter showes good agreement with the recent 
experiments [18]. Although reasonably good agreement 
is found for smaller AO (i.e., larger distance dx)^ the nu- 
merical results strongly deviate from the analytical ones 
for close proximity of the AV to the junction. This faster 
decrease of jc with A^ is due to the currents of the AV, 
which reach and interact with the junction interface. To 
eliminate this effect, we repeated the simulations for a 
sample with twice as large dimensions (i.e., for given A^ 
the vortex is twice further from the junction), which in- 
deed results in a jc very close to the analytic one (see 
filled dots in Fig. [2fb)). Small deviation from the ana- 
lytical results is expected, due to the non-uniform phase 
distribution around the vortex and thus along the junc- 
tion. Very recently Clem [24] studied analytically the 
effect of the nearby vortex on the critical parameters of a 
long planar Josephson junction, who reported several fea- 
tures of the system in common with a recent experiment 

To show properties of the resistive state of the sys- 
tem, and the role of the imprinted phase therein, we 
plotted in Fig. [3] the time evolution of the output volt- 
age V = (p~^ — (p~ for an applied current just above jc, 
for the sample with dx = w/A^ together with snapshots 
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FIG. 4: (color online) (a) I-V curves of the sample of Fig. 
[2] without (solid black curve) and with the pinned Abrikosov 
vortex (symbols) at dx = 500 nm for different displacement 
dy of the pinning site and for two opposite directions of the 
applied current, (b) The critical current jc for two directions 
of the external current as a function of dy (for dx = 500 nm). 



of the distribution of phase and the modulus of the or- 
der parameter. The V{t) characteristic shows periodic 
oscillations (Fig. [3] shows one full period) with voltage 
peaks due to entry /exit of an "antifluxon" in the Joseph- 
son junction (see the insets of Fig. [3]). This Josephson 
antivortex moves in the direction determined by the po- 
larity of applied current, and during its motion forms a 
pair with the pinned AV, which leads to a local minimum 
in the V{t) curve (see inset 3 in Fig. [3j). Note that in 
the absence of the imprinted phase (i.e. pinned AV), the 
dissipation arises from the periodic nucleation and an- 
nihilation of the Josephson vortex- antivortex pairs (not 
shown here) in a similar fashion as in uniform supercon- 
ducting stripes [lH. In our case, the presence of the AV 
does not allow for the formation of a Josephson fluxon- 
antifluxon pair and only the antifluxon contributes to 
the resistive state of our sample, crossing the sample in 
opposite directions for changed polarity of the applied 
current. 

In what follows, we show the dynamics of this anti- 
fluxon in an asymmetrically imprinted phase, obtained 
by shifting the AV off-center in the ^/-direction. Simi- 
larly to Fig. Efa), Fig. Hl^a) shows the I-V curves, but 
now for two different values of the shift dy {dx is flxed). 
In the symmetric case {dy = 0, open and fllled dots in 
Fig. m^a)), jc for both directions of the applied current 
is the same (i.e., = j~) and the output voltage differs 
only by its sign (here we plot \V\) due to the reversal 
of the direction of the antifluxon motion. This situa- 
tion changes entirely for dy ^ 0, i.e. for asymmetric 
distribution of imprinted phase - we observed a smaller 
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FIG. 5: (color online) Vif) characteristics (solid curves) of 
the sample with dx — 500 nm and dy — 500 nm, when biased 
by square-pulsed ac current (dashed curves, right axes) with 
period tp = 500 tql (a) and tp = 100 tql (b). 



needed (negatively polarized) current to induce the an- 
tifluxon motion along the junction (fllled squares in Fig. 
[3Ka)). The reason is the reduced energy barrier for the 
antifluxon entry, which is also reflected as a smaller en- 
try peak in the periodic temporal characteristics of the 
voltage (shown in Supplementary material [26|). On the 
other hand, for the same reason jc increases for positive 
applied current (open squares in Fig. HJ^a)), i.e., jt Jc • 
As a main result, we plotted in Fig. SJ^b) the j+ (fllled 
circles) and j~ (open circles) as a function of the asym- 
metry of imprinted phase (i.e., dy) for dx = 500 nm. The 
difference in jc indicates that for a current amplitude in 
the shaded area of Fig. lU^b) only one direction for mo- 
tion of the antifuxon is possible. In other words, when 
biased by an ac current, our system results in total rec- 
tification of the motion of the Josephson antivortex and 
the corresponding voltage. In our case, the ratchet po- 
tential is created by the imprinted phase fleld (thus not 
by geometry or temporal asymmetry in the junction, as 
is usually the case in vortex ratchets), where the phase 
distribution can be tuned by the position of the AV. For 
that reason, the rectifled voltage, as well as the current 
range for rectiflcation {Aj = j^ —Jc) strongly depend on 
the AV position (see Fig. IHb)). The rectiflcation ceases 
for a symmetric position of the AV {dy = 0), which can 
be restored by a non-zero applied magnetic fleld. 

To further demonstrate the operation of our ratchet 
system, we show its behavior in an ac drive, i.e. a se- 
quence of square current pulses of opposite polarity. Solid 
curve in Fig. [5fa) shows the calculated voltage vs. time 
for the amplitude \j\ = 0.038jo and period tp = SOOtgl 
of the external current. Although the output voltage is 
rather weak (with an average ~ 10 /iF), the ratchet ef- 
fect is clearly seen here. The averaged rectifled voltage, 
which is roughly proportional to the antifluxon veloc- 
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ity, increases with increasing amplitude of the current, 
as also observed in the I-V characteristics of the sample 
in Fig. m^a). Negative voltage at the very moment of 
the current reversal is due to the inertial mass of the an- 
tifiuxon which is already inside the junction. Since the 
voltage rectification in our system is a consequence of 
the directed net motion of the Josephson antifluxon, the 
rectification frequency (i.e., the frequency of the applied 
drive fr = l/tp for which ratchet operates) is restricted 
by the characteristic time scale of the antifiuxon dynam- 
ics - its time needed to pass the entire junction, denoted 
as At. Therefore, if we drive our sample with square 
current pulses with period tp < At, the ratchet behavior 
diminishes and the time averaged voltage becomes zero, 
as shown in Fig. [Sjb) for tp = 100 tgl and At = 145 tgl- 
On a positive side. At decreases with the amplitude of 
applied current, as well as with shortening the junction, 
so that the operating frequency range of our ratchet sys- 
tem increases. In the present calculation, for the chosen 
location of the AV with dx = 500 nm and dy = 500 nm, 
our ratchet rectifies square pulses with frequency up to 
2 GHz, which is much larger than the frequency range of 
ratchets based on moving Abrikosov vortices [3] . We also 
tested our conclusions against sinusoidal current biasing, 
and obtained smaller rectification frequencies compared 
to square pulses with same amplitude, due to difficult 
temporal adjustment between the antifiuxon motion and 
the driving force during the cycle. This fact is already 
well known experimentally (see e.g., Ref. [6]). 

To quantify the efficiency of our ratchet system, and 
compare it to other ratchet realizations, we calculated 
the (long-)time-averaged voltage V^^ for the ac current 
drive relative to time-averaged voltage V^'^ for a dc drive 
of same amplitude, as a function of the frequency of 
the square-pulsed ac current. The ratio r] = y'^^/ydc 
is one of the common measures for the efficiency (see 
Ref. [28] and references therein for other definitions of 
ratchet efficiency), with a maximum of 0.5 for a tempo- 
rally symmetric ac drive. As shown in Fig. EJa), for 
tp > At (for tp < At rectification is not possible), our 
ratchet can have maximal theoretical efficiency which de- 
creases to approximately 0.26 with decreasing frequency 
of the drive. Interestingly, with further decreasing the 
frequency, the maximal efficiency is restored^ and the be- 
havior of efficiency vs. frequency is periodic. This is a 
manifestation of commensurability between the period of 
the drive tp and the time for one Josephson crossing event 
At. As illustrated in Figs. [6fb,c), the maximal efficiency 
is obtained in our ratchet for tp = nAt with n an integer. 
For comparison, earlier Josephson ratchet proposals did 
not have such interesting commensurability effect, and 
also had lower efficiency, e.g. 0.33 in Ref. [27] or 0.22 in 
Ref. [13]. 

Finally, we discuss the efficiency of the phase imprint 
for the realization of a ratchet and compare it to the 
case of an asymmetric potential produced by varying the 
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FIG. 6: (color online) (a) The efficiency 77 (defined in the text) 
of the phase-imprinted ratchet, as a function of the period tp 
of the applied square-pulsed current. (b,c) V{t) characteris- 
tics for periods of the applied drive indicated in (a). 



width of the J J from narrow on one side, to wide on the 
other. The latter geometry should favor a higher volt- 
age for the polarity of the current that drives the fluxon 
from the wide to the narrow end of the junction. How- 
ever, in the absence of any applied magnetic field, we 
observe identical voltage response for both polarities of 
the applied current, i.e. no rectification is possible with- 
out external magnetic field (see Supplementary Material 
for animated data [29]). The reason is the generation 
of either fluxon or antifluxon in the resistive state with 
mirrored spatia l dy namics for the two polarities of the 
applied current [29| . This directly shows that ratchet re- 
alization by imprinted phase is advantageous to earlier, 
geometry assisted vortex ratchets, in the absence of ap- 
plied magnetic fleld. 

In summary, we demonstrated a Josephson vortex ra- 
chet, based on a novel concept - imprinted asymmetric 
phase difference on a Josephson junction. The imprint of 
the phase can be realized by pinning the Abrikosov vortex 
in a suitable location nearby the junction, but even the 
direct phase imprint on the electronic condensate is not 
beyond experimental reach, as demonstrated on Bose- 
Einstein condensates [30] . Our quantum ratchet rectifles 
voltage due to an applied ac current with frequencies up 
to the GHz range (the frequency range is tunable by the 
response time of the sample), even at zero applied mag- 
netic fleld. These properties, and its simple realization, 
make the phase-imprinted Josephson junction an advan- 
tageous fluxon diode in comparison to earlier proposed 
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vortex ratchets. 
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